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We investigated an integral representation of a complex function which was obtained from the 
molecular beam magnetic resonance, and named it as '(/-function' because it was connected with 
Gamma function as a special case. We introduced a physical example of the (/-function from a 
dipole-dipole interaction of rigid polar molecules. 
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Quantum mechanics cannot stand without special functions such as Bessel function or Legendre function. Besides 
these functions, there are so many functions that was originated from mathematics but has been used in physi- 
cal sciences more than mathematics. On the other hand there are still many functions that has been studied by 
mathematician but has been rarely used in other areas. Among them some functions do not have even names. 

In this paper, we will introduce a mathematical function which has not name yet and has been considered as no 
physical application. We call it as (/-function because the gamma function is a very special case of the function. 
Of course it is totally different from the generalized or modified Gamma functions. We'll use small 'g' to avoid the 
confusion from the Meijer's G-function [lj. Then, we will introduce a physical example of the typical case of the 
(/-function, from a familiar classical model of an electromagnetic theory. 
. In various theoretical physics such as a neutron beam velocity distribution function 0, [1[ or a theory of the 
molecular beam magnetic resonance method [H, [H, @, there arise a definite integral function J °° ye~ y ~ z / y dy, where 
z is a complex number. It was the birth of a new complex function containing the Gamma function as a special case. 

In 1937 the (/-function was introduced by Zahn 0] from an absorption coefficients for slow thermal neutrons. The 
fraction of neutrons transmitted is given by 
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g- *W= / ye-'-'l'"'dy. (1) 



Then, the function was extended to m-th order by Laporte [3J . The extended form was given by 
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<Pm(x) = / y m e-y- x ^ 1/2 dy. (2) 
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The order m is restricted to positive integer here. He showed that the integral function is a solution of a third-order 
linear differential equations. Also the convergent expansions were obtained for small values of x. 

Later in 1941 the modern form of the m-th order (/-function was introduced by Torrey Q in a neutron beam velocity 
distribution. The transition probability in radio-frequency spectra contains the function. At result in 1951 the 3rd 
order (/-function was studied mathematically by Kruse and Ramsey [7j . At this time the (/-function was extended into 
complex plane. 

The complete form of the m-th order g-function is defined by 
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9m(z) = / y m e-y 2 - z /ydy, (3) 



where m — 0, 1, 2, and z is a complex number. Interchanging y into 1/x, we obtain another expression of g m {z) as 

poo —l/x 2 —zx 

9m ^=J x m+2 dx - ( 4 ) 

Therefore, the Gamma function is a special case of the (/-function as 

9m (0) = \T (2±1) . (5) 
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Alternatively, 
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5 2m+ i(0) = y. (6) 



The g m has the following three fundamental properties [8j. First, it is a solution of the 3rd order differential equation 

where m = 0, f , 2, ... Second, the derivative reduces single order and change its sign 

= -g m -i(x), (8) 

ax 

where m — 1, 2, 3, .... Third, it satisfies the recurrence relation 

2g m (x) = (m - l)# m _ 2 (a;) + X5 m _ 3 (a;), (9) 

where m = 3, 4, 5, .... More basic properties of the g m in pure mathematical aspects are summarized in the handbook 
by Abramowitz and Stegun Q. 

The g m (z) can be also written as g m (r,9) in the polar coordinate. The complex z is written as z — re~ l6 , where 
r = \z\ and < < 7r/2. Then, the real and imaginary parts are 



Reg m {r, 9) = ^ y m e - y ~ C os — J dy, (10) 

and 

7m 5m (r, 0) = J y m e~y~-— sin ( d y. (11) 

It oscillates very slowly around r axis but converges to zero quickly. 

Furthermore, g m {z), m = 0,1, 2, does not have infinite Taylor series expansion because it is expanded up to the 
rn-th order of z at z = 0. From Eqs. ^ and © the n-th derivatives of g m (z) by z is 

jn / \ poo -l/x 2 -zx 

d -J^l = (_i)« / f dx . (12) 

where n > 1. If we check the value at the origin, 

,. d n g m (z) 

hm = finite constant : n < m. 

|zHo dz™ 

= infinite : n > m. (13) 

Although it has some possibilities of exciting features, there has been no more serious studies in the g m function 
since 1950s. The main reason is that the application area of the function has never been found besides the previous 
known fields. Furthermore, the function was not considered seriously as an extended form of the Gamma function. 
That is, the further studies of the g m was retarded by the restriction of the application area even though it has many 
interesting properties. 

We introduce a classical example of the gi(z) in a complex plane. It comes from the Fourier transformation of the 
time correlation function of dipoles in linear response theory. We begin our story of the function from an example: 
the frequency dependent electric susceptibility of polar molecules. 

We assume a homogeneous isotropic fluid of rigid particles. In general the shape of a rigid molecule is arbitrary, and 
has three different moment of inertia. Also we assume that it has a dipole moment /i(i) in the direction of rotating 
axis. Let u be the applied frequency and r be a relaxation time by collision of the rotating particles. Then, the self 
part of the electric susceptibility has the following form 0, [l(| HH, 03 • 

Xs( " + ; /T) = 1 + i(u> + i/r) r ^^WtA dt . (14) 
Xs{0) Jo M 



(...) denotes an ensemble average over initial conditions or the statistical average in the absence of applied fields. 
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FIG. 1: The radial dependence of gi(r,n/3). 



Let (6i(t),<pi(t),t(ji(t)) be the Euler angles that describe the orientation of particle i at time t with respect to a 
laboratory fixed frame of reference. Also, let and = 1,2,3) be the moment of inertia and angular velocity of 
the three dimensional rigid body to its three principal axes. Then, the angular momentum Li satisfies Li — IiQi. 
The kinetic energy of an arbitrary shape of rotating three dimensional rigid body is composed of translational and 
rotational motion. That is, 
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(15) 



where M is the mass and P is the linear momentum. The f{9, ip) is the angular momentum relation function that is 
defined by 



/(#, ip) = — sin 2 9 sin 2 ip + — sin 2 9 cos 2 ip + cos 2 9 
h h 



(16) 



f(6,ip) varies smoothly in 9 and ip, and it gives us the relation between the rotational kinetic energy and the angular 
momentum. The translational motion of the particle gives just a constant in calculation of the statistical average. 
Therefore, we can neglect P 2 /2M term in further steps. 

For convenience we define the angle difference of the dipole (i at time and t as 7. Then, we can write the 
dipole-dipole correlation function in Eq. (|14[) as 
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cos7(6>, -0, u), 



(17) 



where u is a dimensionless time defined by u = Lt/I^. This conversion makes the cos 7 not include L explicitly 
Applying Eqs. (fTS"|) and (fl7|) to Eq. (fT4)) produces g(9, ip, u) as an essential part of the Fourier transformation. 
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sin9d9 / cos 7 (9, ip, u)g (9, ip,u) dip, 
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where 



g(9^,u)= / L e -W/2/3)£Mi/T-^)W£ di . 
Jo 



(18) 



(19) 



The Z is the partition function and independent of the initial position and time. (3 = l/fc^T, where ks is the 
Boltzmann constant and T is the absolute temperature. 

Eventually, the integral expression of g(9,ip,u) is directly related with the gi(z) in complex plane as 



9i(z) = I ye 
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(20) 
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FIG. 2: The phase dependence of gi(S,9). 

The z is a complex number given by z — z\ — iz-i. Comparing Eq. (| 18[) with Eq. (|19p . the regions of the real number 
z\ and Z2 are decided as Z\ > 0, and zi > 0. We plotted the real and imaginary parts of the gi(r,6). The radial 
dependence at 9 = ir/3 was plotted in FIG. 1, and we see that it oscillates and quickly drops to zero. The phase 
dependence at r = 3 is plotted in FIG. 2, too. 

We named a special function as g m function with order to, and studied some fundamental properties of it. A 
physical example of the <?-function in a complex plane was introduced from classical mechanics and electromagnetic 
theory. 

Although, some functions has been regarded as not useful application areas, it is not actually true. We have not 
found the application area yet. It is interesting enough to study connecting with the Gamma function. Merging to 
of the radial part is quicker than Bcssel functions and Neumann functions, and that may applicable to a new integral 
transformation. We hope that it covers various range of physics that previous ordinary Gamma function could not. 
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